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A variable combination of realistic and random two-body interactions allows the study of collective 
properties, such as the energy spectra and B(E2) transition strengths in 44 Ti, 48 Cr and 24 Mg. It 
is found that the average energies of the yrast band states maintain the ordering for any degree of 
randomness, but the B(E2) values lose their quadrupole collectivity when randomness dominates 
the Hamiltonian. The high probability of the yrast band to be ordered in the presence of pure 
random forces exhibits the strong correlations between the different members of the band. 

PACS numbers: 21.10.Re, 05.30.Fk, 21.60.Cs, 23.20.-g 
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I. INTRODUCTION 

Nuclei display very regular spectral patterns. Low en- 
ergy states in medium- and heavy-mass even-even nuclei 
allow their classification in terms of seniority, anharmonic 
vibrator and rotor nuclei [jj] , according to the ratio of the 
excitation energies of the states 4i and 2i. While this 
regular behavior has been usually related with specific 
forces, the investigation of the energy spectra with ran- 
dom interactions y, |3| has shown that many body states 
are strongly correlated even in the presence of random 
two-body interactions. Random interactions in bosonic 
Hilbert spaces, like those used in the IBM and the vi- 
bron model, exhibit a large predominance of vibrational 
and rotational spectra, strongly suggesting that in boson 
spaces collectivity is an intrinsic property of the space of 
nuclear states ^, 0, ||] . Recently several studies were per- 
formed with random (tbre) and displaced random ensem- 
bles in order to simulate realistic systems |?], ||, in par- 
ticular to investigate the dominance of 0+ states 

In the nuclear shell model the transition between ran- 
dom and collective behavior in the energy spectra of 
20 Ne generated by two-body forces was addressed in [fl2| . 
Collectivity was generated with a quadrupole-quadrupolc 
force, while a residual random interaction was included 
in the Hamiltonian in order to study its consequences on 
the system's spectroscopic properties. Both the eigen- 
value distribution and the overlap between the SU(3) and 
calculated wave functions exhibit the smooth path from 
a Hamiltonian dominated by the collective force to a ran- 
dom, non-collective, one. 

The present work aims to extend these ideas by study- 
ing the transition from a realistic parametrization of the 
two-body force to a purely random one, in complex sys- 
tems like 24 Mg, 44 Ti and 48 Cr. The realistic interactions 
we have chosen are the universal Wildenthal jl3| inter- 
action for the sd-shell, and the KB3 [Q for the fp-shell. 
The probability for each state in the yrast band to follow 
a sequence where the higher energies correspond to the 



states with the larger angular momentum is studied by 
varying the mixing between realistic and random forces 
in the Hamiltonian. The evolution of the average ener- 
gies for each member of the band, as well as its B(E2) 
values, is also reported. 



II. BAND STRUCTURE 



The combination of random and realistic interactions 



is taken as |12 

H = aH c + bH R , with a + b = 1, 



(1) 



where He is a realistic Hamiltonian and Hr is a two 
body random ensemble. Both He and Hr are written as 



H 



E 



V JT 7 \ 
v rstu Zj rs,JT 



->tu,JT 



(2) 



in terms of scalar products of the normalized pair 
creation operators Zj. JT and its Hermitian conjugate 
Zm,jt, where r, s, . . . specify sub-shells associated with 
individual orbits, and J, T the coupled angular momen- 
tum and isospin. For the realistic Hamiltonian V r J s J u is 
the Wildenthal or K B3 interaction. For the random case 
the V r J s J u matrix elements are taken from a two body ran- 
dom ensemble (TBRE), i.e., to be real and normally dis- 
tributed with mean zero and width a for the off-diagonals 
and v2c for the diagonals. The values of the width a are 
taken from the realistic interactions: 1.34 MeV and 0.60 
MeV, for the sd- and fp-shells, respectively. The param- 
eter b is varied from to 1, to cover the different mixing 
from the realistic interaction to a pure random force. 



A. The fp shell 

As a first example we take 48 Cr, a rotational nucleus 
which has been widely studied in full shell model cal- 
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culations. In the present case the diagonalization is per- 
formed in the shell with the code Antoine |16|] , using 
a KB3 interaction jllj] without single particle energies. 
Fig. 1 shows the average energies Ej of the lowest en- 
ergy state for each angular momentum J= 0, 2, 16 h, 
the yrast band, calculated for 960 samples of the random 
interaction. 
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FIG. 1: Average energy of 960 runs for each angular momen- 
tum of the ground band in 48 Cr. 

The vertical axis displays the average energy of each 
state, as a function of the mixing parameter b. At the 
left hand side the realistic energies are shown, with its 
distinctive rotor pattern. As b increases to the right, 
the order between the different members of the bands is 
maintained, but their relative separation changes. The 
average energies in the right hand side, the pure random 
Hamiltonian, still exhibit a band structure but have lost 
their quadrupole collectivity, as discussed below in con- 
nection with their B(E2) values. 

The evolution of the average ground state band of 44 Ti, 
calculated in the full pf-shcll, is shown in Fig. 2 as a func- 
tion of the mixing parameter b. It can be seen that the 
relative ordering of the states with different angular mo- 
mentum is maintained, but their relative separations vary 
significantly in the transition from the realistic to the 
random Hamiltonian. However, up to b = 0.6 the change 
is mostly a scale variation, with the nearly equidistant 
structure of the band keeping its form. The energy spec- 
tra evolves from a vibrational equidistant form at b = 
to a mixed, yet ordered, spectrum for pure random forces. 
The average ground state energy increases with the mix- 
ing to a maximum value, with a slight decrease for 6=1. 

Similar patterns of evolution of the energy centroids 
for each angular momentum, as function of the mixing 
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FIG. 2: Average energy of 960 runs for each angular momen- 
tum of the ground band in 44 Ti. 



parameter, were found in several other nuclei like 46,48 Ca 
and 46 Ti. While for 46 Ti, the centroid of the 3=2 state 
is lower in energy than the J=0 state for a pure random 
Hamiltonian, in general for these nuclei and others in 
the sd-shell the average energy ordering is conserved, and 
there is a gradual change in the energy spectra. 

The width a associated with the lowest average energy 
for each angular momentum, calculated as the square 
root of the variance, is shown in Fig. 3 for the gs-band of 
48 Cr, as a function of the angular momentum, for differ- 
ent values of the mixing parameter b. The energy width 
increases with the mixing, from 1.5 MeV for b — 0.2 to 
7-9 MeV for b = 1.0. Given that, as b increases, these 
widths become larger than the gaps between the average 
energies, they indicate a strong mixing between different 
bands, an effect closely associated with the lack of col- 
lectivity discussed bellow. As a function of the angular 
momentum all the widths are essentially flat, showing a 
moderate increase for the largest angular momenta. 

The energy widths are larger than the relative energy 
differences between adjacent states in the g.s. band, im- 
plying that there are individual cases in which the spec- 
trum is not ordered. However, the number of runs with a 
completely ordered spectrum is large even for pure ran- 
dom interactions, as is discussed in detail bellow. 

B. The sd-shell 

Calculations in the sd-shell can be done for the full 
shell and an arbitrary number of active particles. The 
realistic interaction used in this case is the universal 
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FIG. 3: Average energy width of 960 runs for each angular 
momentum of the ground band in 48 Cr. 



FIG. 4: Average energy of 960 runs for each angular momen- 
tum of three bands in 24 Mg, see text for labels description. 



Wildenthal interaction (USD pfl). The nucleus 24 Mg 
offers a rich enough system, where three bands can be 
studied simultaneously. They correspond, for b = 0, to 
the ground state (gs), (3- and 7- bands. The first two 
start with J=0 and contain only even-J states, while the 
7- band starts with J = 2 and includes states with both 
even and odd angular momenta, up to J = 8. Following 
the evolution of the average energy levels as the mixing 
parameter b increases, we reconstruct the three bands 
for each value of b. Nevertheless, for a pure random in- 
teraction one can show that the states have lost their 
quadrupole collectivity (see bellow), and for this reason 
they do not form bands in the usual sense. 

In Fig. 4 the average energies for each angular mo- 
mentum, for the ground, (3- and 7- bands in 24 Mg, are 
presented for different values of the mixing parameter b. 

The most relevant features observed in Fig. 4 are that 
the three bands evolve in a similar way, keeping their in- 
ternal order as well as the relative separation between the 
bands. The correlation between increasing energy and 
angular momentum is strictly followed in all the aver- 
aged bands. The absolute ground state energy becomes 
larger as the mixing parameter increases, following the 
same pattern found in the fp-shell. 



the highest the energy for a given band. To do so we 
counted the number of runs where the first J=0 state is 
the ground state, or the first excited state in the band, 
the second excited state, and so on, in the 960 runs. We 
did the same for the states with angular momenta J = 2, 
4, ... In Table I the percentages for 48 Cr with a mixing of 
b = 0.4 are listed. All the states have a probability of at 
least 90% to occupy its physically expected place, while 
the dispersion is very small. 



J+ 


1st 


2nd 


3rd 


4th 


5th 


6th 


7th 


8th 


9th 





96.56 


1.04 


1.25 


0.52 


0.42 


0.10 


0.00 


0.10 


0.00 


2 


1.56 


91.88 


5.31 


0.62 


0.21 


0.21 


0.10 


0.00 


0.10 


4 


0.21 


5.62 


92.29 


1.04 


0.42 


0.21 


0.21 


0.00 


0.00 


6 


0.52 


0.94 


0.42 


92.50 


4.27 


0.94 


0.21 


0.00 


0.21 


8 


0.73 


0.21 


0.21 


3.75 


93.02 


1.04 


0.31 


0.52 


0.21 


10 


0.00 


0.00 


0.10 


0.83 


0.52 


95.21 


2.29 


0.83 


0.21 


12 


0.00 


0.00 


0.21 


0.62 


0.83 


2.08 


96.04 


0.21 


0.00 


14 


0.21 


0.21 


0.10 


0.10 


0.31 


0.10 


0.52 


97.40 


1.04 


16 


0.21 


0.10 


0.10 


0.00 


0.00 


0.10 


0.31 


0.94 


98.23 



TABLE I: Probability for each angular momentum state to be 
in the correct position in the ground state band with b — 0.4 
for 48 Cr. 



III. ORDERING PROBABILITIES 

In order to further study these systems, we have ana- 
lyzed the probability that each state in the band has the 
usual ordering i.e, the largest the angular momentum, 



Table II lists the probability for states of each angular 
momentum in 48 Cr to occupy the indicated place for a 
purely random Hamiltonian, i.e. 6=1. The probabilities 
of being in its expected place run from 44% for J = 2 
to 83% for J = 16. States with J = 0, 2 and 4 are 
those which more often fail to occupy their place, as they 
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tend to exchange positions. These results represent an 
extension of previous studies concerning the probability 
of each state to be the ground state, listed in the first 
column U). 



J+ 


1st 


2nd 


3rd 


4th 


5th 


6th 


7th 


8th 


9th 





60.31 


8.85 


9.27 


4.48 


4.90 


2.50 


1.98 


2.29 


5.42 


2 


17.40 


43.54 


16.98 


4.17 


3.96 


2.60 


3.12 


5.00 


3.23 


4 


6.56 


27.71 


47.08 


6.46 


2.40 


2.71 


4.38 


2.19 


0.52 


6 


2.40 


5.42 


9.79 


64.58 


6.25 


5.31 


3.12 


1.77 


1.35 


8 


3.12 


3.65 


4.69 


9.06 


69.58 


6.46 


1.15 


1.56 


0.73 


10 


1.35 


2.50 


3.44 


5.94 


7.19 


74.69 


3.54 


0.62 


0.73 


12 


1.98 


1.88 


5.52 


2.50 


3.23 


3.65 


79.06 


1.98 


0.21 


14 


1.15 


5.62 


1.98 


1.67 


0.94 


1.25 


2.40 


80.10 


4.90 


16 


5.73 


0.83 


1.25 


1.15 


1.56 


0.83 


1.25 


4.48 


82.92 



TABLE II: Probability for each angular momentum state to 
be in different positions, for the 48 Cr ground state band with 
6=1.0. 



Table III displays the probability for each angular mo- 
mentum state to be in a given position, in the 44 Ti ground 
state band with b = 1.0. These probabilities have values 
between 44% for J=2 to 95% for J = 12. The fact that 
the states with J=2 and 4 have less than 50% probabil- 
ity of occupying in their places is strongly connected with 
the closeness of their average energies, shown in Fig. 2 
for 6 = 1. 

J+ 1st 2nd 3rd 4th 5th 6th 7th 
46.25 19.27 20.42 6.98 5.62 1.35 0.10 
2 16.56 39.90 24.38 11.56 5.83 1.25 0.52 
4 21.04 25.83 40.10 9.38 3.02 0.21 0.42 
6 8.44 10.42 10.10 64.58 5.42 1.04 0.00 
8 6.56 3.96 4.27 5.83 76.77 2.08 0.52 
10 0.10 0.52 0.52 1.25 2.60 92.08 2.92 
12 1.04 0.10 0.21 0.42 0.73 1.98 95.52 

TABLE III: Probability for each angular momentum to be 
in different positions, for the 44 Ti ground state band, with 
b = 1.0. 



Tables IV, V and VI display the probability that states 
with different angular momenta in 24 Mg have to occupy a 
given place in each band for b = 1.0, for the ground state, 
/3- and 7- band, respectively. In most cases the diagonal 
probability, i.e. the probability that each state occupies 
its expected place, is larger than 50%. The exceptions 
are the states with J = 2 and 3 in the 7-band, and those 
with J= and 2 in the /3-band, whose probabilities lie 
between 38% and 48%. 



J+ 1st 2nd 3rd 4th 5th 

56.56 12.40 12.29 6.77 11.98 

2 18.12 52.81 13.12 9.79 6.15 

4 9.06 19.58 63.85 5.52 1.98 

6 4.38 11.25 6.98 71.25 6.15 

8 11.88 3.96 3.75 6.67 73.75 

TABLE IV: Probability for states belonging to the ground 
state band in 24 Mg, with b = 1.0, to occupy different posi- 
tions. 



J+ 1st 2nd 3rd 4th 5th 6th 7th 

2 45.21 22.92 7.50 5.94 4.79 6.46 7.19 

3 27.40 38.33 13.54 6.15 4.17 6.04 4.38 

4 7.40 14.90 50.62 15.94 7.29 3.02 0.83 

5 7.71 12.08 15.10 58.85 5.00 1.04 0.21 

6 3.44 3.75 6.25 7.60 66.25 11.88 0.83 

7 5.10 4.27 4.17 3.96 10.83 64.38 7.29 

8 3.75 3.75 2.81 1.56 1.67 7.19 79.27 

TABLE V: The same as in Table 4, for the 7-band. 



IV. COLLECTIVITY 

A sensitive measure of the quadrupole collectivity of 
the system is the B(E2,2i — * gs ) transition strength. 
The distribution of B(E2) strengths in 48 Cr, for four dif- 
ferent values of the mixing parameter b, is shown in Fig. 
5. Its experimental value is 230 e 2 b 2 . For b = 0.2, shown 
in insert a), the distribution is concentrated around the 
experimental B(E2) value. For b = 0.4, insert b), the 
increase in the random components of the Hamiltonian 
leads to an important fragmentation of the B(E2) in- 
tensity with four clusters, one near zero, the second one 
around 75 e 2 6 2 , the next one near 170 e 2 6 2 and the last 
one close to the measured value. For b = 0.6, insert c), 
most of the B(E2) values are very small, with some in- 
tensity at the collective B(E2) values. Finally, in insert 
d) the distribution of B(E2) values for a purely random 
Hamiltonian is shown. It is strongly concentrated at very 
small values, showing a complete lack of quadrupole col- 
lectivity, in consonance with the findings of Ref . |l^| . 

Fig. 6 displays the distribution of B(E2,2 — > 0) values 
for 44 Ti. Insert a) shows the results for b = 0.2, with 
a narrow distribution around the measured value of 147 
e 2 b 2 . In insert b) the distribution for b — 0.4 is presented, 
which is concentrated around the same value. For b = 
0.6, Fig. 5 c), the distribution is still concentrated around 
the collective B(E2) values despite the dominance of the 
random component in the Hamiltonian. However, this 
collectivity is completely lost when pure random forces 
are employed, as shown in Fig. 5 d). 
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J+ 1st 2nd 3rd 4th 5th 

46.67 20.73 13.54 8.44 10.62 

2 28.02 47.40 12.81 8.44 3.33 

4 13.23 19.27 65.62 1.88 0.00 

6 6.35 9.38 5.94 77.50 0.00 

8 5.73 3.23 2.08 3.75 85.21 

TABLE VI: The same as in Table 4, for the /3-band. 
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FIG. 5: Probability densities for B(E2,2 -> 0) for 48 Cr. 

A useful indicator of collectivity is the energy ratio 
E 4t 



R = 



E 2+ 



(3) 



The value R = 2 is associated with a harmonic oscilla- 
tor spectrum, while R = 3.3 characterizes a rigid rotor 
structure. 

Fig. 7 shows the distribution of energy ratios for 48 Cr. 
The case b = 0.2, shown in insert a), does exhibit the 
actual rotor behavior of this nucleus. This feature re- 
mains dominant for b = 0.4, insert b), while for b = 0.6 
the distribution of energy ratios is wide and peaked at 
R=l. For b = 1.0, shown in Fig 7 d), the distribution 
is very wide, with a clear dominance of the R=l ratio, 
in correspondence with the near degeneracy of the states 
with J= 2 and 4 for b = 1, as shown in Fig. 1. 

4 Ti has a structure closer to a harmonic oscillator 
than to a rotor. This feature can be seen in Fig. 8 
a), which displays the distribution of energy ratios for 
b = 0.2, narrowly concentrated around R = 2. For b = 0.4 
and 0.6, inserts b) and c), the vibrational structure is 
wider but well defined. For a pure random interaction, 
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FIG. 6: Probability density for B(E2,2 -» 0)for 44 Ti. 
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FIG. 7: Probability densities for the energy ratios R for 48 Cr. 



Fig. 8 d), the distribution is peeked at R = 1, reflecting 
the near degeneracy of the average energies Ei and £4. 
The displacement of the most probable energy ratio R to 
1 is accompanied by a lack of quadrupole coherence, in 
consistent with the previous analysis of the B(E2) tran- 
sition strengths. 
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FIG. 8: Probability densities for energy ratios R for Ti. 



FIG. 9: Percentage of ordered bands vs b for Cr,Ti and Mg. 



V. CORRELATION AND COHERENCE. 

Having discussed the evolution of the energy centroids 
as a function of the mixing parameter b, and the probabil- 
ity for each state to occupy a certain position in the band, 
it is natural to study the probability that the ground 
state band has all its states properly ordered. This prob- 
ability, given as the percentage of results from the 960 
runs which are properly ordered, is shown in Fig. 9 for 
48 Cr, 44 Ti and 24 Mg. In the three cases it is apparent 
that the spectrum is always ordered when the realistic in- 
teraction dominates the Hamiltonian, and that the prob- 
ability for the ground state band to be ordered decreases 
to about 35% for a fully random Hamiltonian. 

It is worth emphasizing that if these probabilities were 
independent of each other, the probability to find a com- 
pletely order band would be the product of the proba- 
bilities for each state, i.e., the product of the diagonal 
elements of the matrices shown in Tables II, III and IV. 
However, these products are all smaller than 0.04, far 
smaller than the probabilities of 25 - 45 % found in these 
calculations. These results show that the many body 
states obtained with two-body random forces are strongly 
correlated. This correlation in their energies is not, how- 
ever, associated with quadrupole coherence, which is lost 
for pure random forces fl2|| . The subtle connection be- 
tween correlation and coherence requires further analysis. 

The correlation r for each value of the mixing param- 
eter &, between the states with energies Ej and Ej+2 in 
the sample, can be quantitatively calculated as |19| 



cov(Ej,E J+2 ) 



where 



cov(Ej,E, 



J+2) 



J2i ( E Ji - E j)( E (J+2) i - E J+2 ) 



(•5) 



Here, cov(Ej,Ej + 2) is the covariance of the two distri- 
butions, a is the square root of the variance cov{Ej, Ej), 
and n is the size of the space. If r = 1 the distributions 
are fully correlated, while if r « there is no correlation 
between them. In table VII the correlation r for pairs of 
energies in the ground state band of Cr are listed . 

J, J + 2 0.2 0.4 0.6 0.8 1.0 



0+,2+ 0.966 0.928 
2+,4+ 0.996 0.991 
4+,6+ 0.992 0.977 
6+,8+ 0.991 0.989 
8+, 10+ 0.987 0.989 
10+, 12+ 0.985 0.982 
12+, 14+ 0.992 0.989 
14+, 16+ 0.992 0.992 



0.937 0.949 0.991 

0.991 0.991 0.986 

0.974 0.981 0.980 

0.992 0.993 0.993 

0.989 0.989 0.987 

0.982 0.984 0.985 

0.987 0.986 0.986 

0.990 0.989 0.988 



O-EjO-E., 



(4) 



TABLE VII: Correlations r between neighbor states in 48Cr. 



The presence of strong correlations between the dif- 
ferent wave functions obtained with two-body random 
forces, suggests that the many-body states could be well 
approximated by a small number of configurations which 
may correspond to definite shapes, as was found for 
bosonic models. This would imply that the very large 
number of shell-model many-body states would be lim- 
ited or constrained by the geometry imposed by the exis- 
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tence of a two-body Hamiltonian, even for the case that 
its components are randomly selected |llf . 

VI. THE DTBRE CASE 

In || it was shown that a displaced two body random 
ensemble (DTBRE) gives rise to coherent rotor patterns. 
To complement the present study we analyze in this sec- 
tion the transition from the realistic KB3 interaction to a 
DTBRE in the ground state band of 48 Cr. The DTBRE 
correspond to matrix elements with a normal distribution 
centered at c = —1.0 and width a = 0.6. Fig. 10 shows 
the evolution of the average energies in the ground-state 
band of 48 Cr as a function of the mixing parameter b. 
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FIG. 10: Average energy of 960 runs for each angular mo- 
mentum of the ground band in 48 Cr with a DFBRE mixing. 

At variance from the results shown in Fig. 1, 2 and 
4, when the random component of the interaction in- 
creases the absolute energies continuously decrease. As 
expected, the rotor structure survives the transition from 
realistic to random interactions with nearly no changes. 

In Table VIII we show the probability for states be- 
longing to the ground state band in 48 Cr, to occupy its 
expected place. Notice that even for b = 1.0, when the 
pure displaced random ensemble is employed, the proba- 
bility for each state to be in its place is larger than 88%, 
thus exhibiting a high ordering. The percentage of fully 
ordered bands for each mixing is equally large: 100% 



for b = 0.2, 89% for b = 0.6 and 70% for b = 1.0. If 
the centroid of the DTBRE is displaced to more negative 
values, like c = —3.0, 100 % of the bands are ordered 
||. The distribution of B(E2) probabilities exhibits a 
clear presence of quadrupole coherence for c = — 1 . They 
are concentrated in a narrow peak around the collective 
B(E2) values for c= -3 §. 



J+ b = 0.2 


0.4 


0.6 


0.8 


1.0 





99.58 


99.79 


97.71 


92.08 


88.12 


2 


98.75 


99.27 


95.00 


86.04 


80.31 


4 


98.96 


99.38 


95.21 


86.46 


79.90 


6 


99.27 


99.58 


96.35 


91.46 


87.50 


8 


99.27 


99.69 


98.02 


94.69 


91.15 


10 


99.58 


99.90 


98.96 


96.77 


94.38 


12 


99.69 


99.90 


98.65 


96.77 


95.00 


14 


99.58 


99.58 


94.79 


90.62 


87.92 


16 


99.58 


99.69 


95.00 


90.42 


88.33 



TABLE VIII: Probabilities for states belonging to the ground 
state band in 48 Cr, to occupy their ordered place, for a DT- 
BRE mixing. 



VII. SUMMARY AND CONCLUSIONS 

The average energies of states with different angular 
momentum preserve their ordering inside the band when 
the Hamiltonian is changed smoothly from a realistic to a 
random one. Ground state energies increase as a function 
of the mixing parameter. The quadrupole collectivity is 
lost when the Hamiltonian is dominated by random two- 
body forces, and the probability that the ground state 
band remains ordered diminishes to 25-45% in the ran- 
dom limit, which is anyway far larger than the product 
of the probabilities for each state to be in its place, thus 
exhibiting the strong correlations between the different 
wave functions. On the other hand, when displaced two- 
body random ensembles are employed, the average en- 
ergies decrease as the random component increases, and 
the rotor pattern remains unchanged. 
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